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Abstract 



The stochastic cellular automaton of Rule 18 defined by Wolfram [Rev. 
Mod. Phys. 55, 601 (1983)] has been investigated by the enhanced coherent 
anomaly method. Reliable estimate was found for the (3 critical exponent, 
based on moderate sized (n < 7) clusters. 
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Calculating critical exponents of second order phase transitions is a challenging problem. 
For nonequlibrium systems, generalization of equilibrium statistical physics methods is under 
developement. Among the most notable analytical tools are series expansion [Q], transfer 
matrix diagonalization || and the mean- field renormalization-group method 0. 

In a series of earlier papers we have shown how the generalization of the mean- field 
technique with appropriate extrapolation can be used to describe the critical properties of 
cellular automata (CA) phase transitions. 

The generalized mean-field approximation (GMF) first proposed for dynamical systems 
by Gutowitcz et al. |7J and Dickman || is shown to converge slowly at criticality. In this 
method we set up equations for the steady state of the system based on n-point block 
probabilities. Correlations with a range greater than n are neglected. By increasing n from 
1 (traditional mean-field) step by step we take into account more and more correlations 
and get better approximations. The GMF approximation can be used as the basis of the 
coherent anomaly method (CAM) calculation, and it gives a reasonably good (3 exponent for 
a dynamical system with large (n > 10) ||. In this letter I show how an improved version 
of the CAM proposed very recently [K| works on cellular automata. 

The essence of the CAM is that the solution for singular quantities at a given (n) level 
of approximation {Q n (p)) in the vicinity of the critical point is the product of the classical 
singular behavior multiplied by an anomaly factor (a(n)), which becomes anomalously large 
as n — > oo (and — » p c ): 

Q n ~a(n)(p/p n c -ir- (1) 

where p is the control parameter and uj c i is the classical critical index. The divergence of 
this anomaly factor scales as 

a(n) ~ (p n c - p c r^ (2) 

thereby permitting the estimation of the true critical exponent u, given a set of GMF 
approximation solutions. However such an estimation depends to some extent on the choice 
of the independent parameter (p <-> 1/p). To avoid this a corrected CAM was proposed fLOf , 
based on a new parameter: 

Sn = (Pc/P n c ) 1/2 - (P n c /Pc) 1/2 (3) 

such that Eq. (Q) is invariant under p p~ x . This parametrization gives better estimates 
for the critical exponents of the 3-dimensional Ising model |TD . 



My target system for this kind of calculation was the one-dimensional, stochastic Rule 
18 CA [|12|| . This range-one cellular automaton rule generates a 'V at time t only when the 
right or the left neighbor was 'V at t — 1 : 

t-1 : 100 001 

t : 1 1 



with probability p. In any other case the cell becomes '0' at time t. The order parameter is 
the concentration (c) of 'l'-s. For p < p c the system evolves to an absorbing state (c = 0). 
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For p > p c a finite concentration steady state appears with a continuous phase transition. 
This transition is known to belong to the universality class of directed percolation (DP) [fl3 



At t — > oo the steady state can be built up from '00' and '01' blocks [|14]]. This permits one 
to map it onto stochastic Rule 6/16 CA with the new variables '01' — 1' and '00' — >' 0': 



t-1 : 00 01 10 11 

t : 1 10 

and the GMF equations can be set up by means of pair variables. In an earlier work [|J this 
was performed up to the order n = 6, and Pade extrapolation was applied to the results. 
Our best estimate for critical data was p c = 0.7986 and (3 = 0.29. 

Now, I have extended the GMF calculations up to n = 7 (see Table |) with the help of 
the symbolic Mathematica software. This required the setting up and solution of a set of 
nonlinear equations of 72 variables. I obtained p 7 c = 0.7729, which is still 5% off the result 
obtained by steady state simulation: p c = 0.8086(2) |15| or from the more accurate time 
dependent simulation data: p c = 0.8094(2) ]TB[ . 

The CAM analysis of (a(n),S n ) data was done, taking into account the correction term 



a\n) 



bSt^ + cSt^ 1 (4) 



and examining the stability of the solution by omitting different points from the (n = 1, .., 7) 
data set. For the fitting (3 c i = 1 and p c = 0.8094 were used. As was pointed out in Ref. 
|TT| the CAM data may contain departures from ideal scaling; moreover there is no clear 



dependence on the order of the approximations. I found relatively stable estimates using 
the correction formula (Q) on the data set with the omission of the n = 3 point. The n = 3 
approximation result does not fit into the (log(5 n ) — log(a(n)) curve either (see Fig. ). Table 
JH shows the stability of the results, with the mean (3 = 0.2796(2) calculated from them. 
This compares very well with the value f3 = 0.2769(2) obtained by Dickman and Jensen |I[ 
from series expansion. If the CAM calculation based on p or 1/p independent variables the 
results differ by ±0.005 from the present enhanced version (3 estimates. 

Another critical model with non-DP universality, the nonequilibrium kinetic Ising model, 
has been examined with the enhanced CAM method, and the (3 exponent estimate is in 



agreement with the simulation results [17 



The conclusion of this study is that the enhanced version CAM method with careful 
data analysis gives good estimates for the critical exponent for moderate n < 10 level 
GMF approximations. Calculation of the n = 5, 6, 7... GMF approximations is possible on 
moderate sized workstations. The solution of the n = 7 level approximation took about 10 
hours CPU time on a SUN Sparc- 10. This provides an efficient analytical tool for exploring 
universalities of nonequlibrium systems. 

This research was partially supported by the Hungarian National Research Fund (OTKA) 
under grant numbers T-4012 and F-7240. 
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FIGURES 

Result obtaining by applying the improved CAM method on n-pair (n = 1, .., 7) approx- 
imation data. The logarithm of the anomaly coefficient a(n) is plotted versus the logarithm 
of the improved independent variable 5 n . Fitting was done according to Eq. (|J). 
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TABLES 

TABLE I. GMF calculation results for pair approximation data 



n a(n) 

1 0.5000 0.5000 

2 0.6666 1.5000 

3 0.7094 2.3484 

4 0.7413 2.8816 

5 0.7543 3.5345 

6 0.7656 4.2545 

7 0.7729 4.8463 



TABLE II. CAM calculation results for pair approximation data 



data 


P 


1-2-4-5 


0.273 


1-2-4-5-6 


0.271 


1-2-4-5-6-7 


0.282 


1-2-4-5-7 


0.285 


1-2-4-6-7 


0.275 


1-2-5-6-7 


0.310 


1-4-5-6-7 


0.275 


2-4-5-6-7 


0.266 


mean 


0.2796(2) 



Pade extrapolation, ref. M 0.29 



simulation, ref. g§ 0.285(5) 



series expansion for DP ref. ffl 0.2769(2) 
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